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Abstract 

We consider the nonlinear stability of a fully three-dimensional boundary layer 
flow in an incompressible fluid and derive an equation governing the nonlinear de- 
velopment of a stationary cross-flow vortex The amplitude equation is a novel 
mtegro-differential equation which has spatial derivatives of the amplitude occur- 
ring m the kernal function. It is shown that the evolution of the cross-flow vortex is 
strongly coupled to the properties of an unsteady wall layer which is in fact driven 
by an unknown slip velocity, proportional to the amplitude of the cross-flow vor- 
tex. The work is extended to obtain the corresponding equation for rotating disk 
flow. A number of special cases are examined and the numerical solution for one 
of cases, and further analysis, demonstrates the existence of finite-distance as well 
as focussing type singularities. The numerical solutions also indicate the presence 
of a new type of nonlinear wave solution for a certain set of parameter values. 
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1 Introduction 


One of the earliest experimental and theoretical investigations of the stability of thiee- 
dimensional (3D) boundary layers was conducted by Giegory. Stuart & Walker (1955). 
(hereafter referred to as GSW). The boundary layer flows studied were the flow over a 
rotating disk, and the flow over swept wings. Using a china-clay visualisation technique 
they were able to demonstrate the presence of a highly legular, stationary, pattern of 
vortices spaced equally around the disk, or along the surface of the wing In addition, 
with the aid of a microphone probe, they were able to detect travelling waves close 
to the surface of the disk. Stuart in GSW suggested that the instabilities could be 
explained in terms of the inflexional character of the effective mean velocity profile in 
certain directions, ( the term effective mean velocity profile refers to a certain linear 
combination of the streamwise and spanwise velocity components) According to his 
suggestion the stationary pattern was that associated with the inviscid instability of the 
velocity profile which had a zero at a point of inflexion. The non-stationary, or travelling 
wave pattern, could also be explained in terms of the inviscid instabilty of the mean flow. 
These instabilities are now commonly referred to as cross-flow instabilities and their 
importance stems from the fact that they arise naturally in many fully three-dimensional 
boundary layer flows of practical importance, such as in the flow past swept aircraft 
wings, rotating flows, and so on. 

The experimental observations of GSW have been confirmed in many subsequent 
experiments on rotating disk flow, see for instance Malik et al. (1981), Kohama (1984). 
Kobayashi et al. (1980), Kohama et al. (1991), as well in experiments on swept wings 
and cylinders, Poll (1985), Michel et al. (1985), Mullei Bippes (1988), Saric el. (1989). 
An extensive review of many aspects of the instability of 3D boundary layers including 
cross-flow instability, and covering many of the early investigations, may be found in 
Reed & Saric (1989), see also Arnal (1986). 

The more recent experiments on cross-flow instability have considered nonlinear effects 
and in Kohama (1984), Kohama et al (1991) an explosive secondary instability, which 
takes the form of ring-like vortices wrapped around the primary cross-flow vortex, has 
been observed. In the experiments of Muller & Bippes (1988) it has been found that 
the stationary as well as the non-stationaxy vortices attained a nonlinear saturation 
amplitude. A similar phenomenon has been observed in the direct numerical simulations 
of Meyer & Kleiser (1988), and also Malik (1986), Malik Li (1992), Malik et. al. (1994). 
using more approximate methods based on the parabolised stability equations approach 

There are very few self-consistent theoretical investigations of the stability of cross- 
flow vortices. Hall (1986) extended Stuart’s analysis to compute correction terms for 
the number of vortices as well as the inclination of the stationary vortices in the limit 
of large Reynolds numbers. He was able to explain the existence of another branch 
of the neutral curve, calculated previously by Malik (19S6). as being associated with a 
mode with zero shear stress. Mackerrel (1987) extended Hall's work to derive a weakly 
nonlinear amplitude equation for the zero shear stress mode. Bassom & Gajjar (1989) 
studied the linear and nonlinear neutral stability properties of long wavelength non- 
stationary cross-flow vortices using scalings appropriate to the upper-branch stability of 
a two-dimensional boundary layer. Recently Gajjai (1994). (1995), has extended the 
Bassom &; Gajjar (1989) work to compressible flows and looked at non-neutral modes 

The mam objective of this paper is to present a self-consistent theoretical description 
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of the nonhneai evolution of a stationary cross-flow voitex. In companion papers Gajjar 
& Arebi (1995) this vvoik is extended to the non-stationai y case and in Gajjar k Sibanda 
(1995) to compiessible flows We study a general 3D boundaiy layer flow and consider 
a flow direction in which the effective velocity profile has a zero at a point of inflexion, 
giving rise to the stationary cioss-fiow voitex. and include modulation m time and space. 
The stability analysis leads to a novel mtegro-differential equation for the amplitude of 
the cross-flow vortex with spatial derivatives of the amplitude occurring in the kernal 
function. 

The starting point for the analysis is to consider a stationary cross-flow vortex in a 
neutral or near neutral state Thus with small growth rates, and since the basic instabilty 
mechanism is Rayleigh instability, GSW (1955), we may appeal to the ideas of unsteady- 
nonlinear critical layer theory used successfully in many studies of the instabilities of 
planar shear and boundary layer flows, see for example Hickernell (1984), Goldstein k 
Leib (1989), Leib(1991), Goldstein k Choi (1989), Wu(1992), Wu (1994). The application 
of unsteady nonlinear critical layer theory to shear flow and boundary layer instabilities, 
and a discussion of the properties of some of the mtegro-differential amplitude equations 
which arise in this type of work, may be found in the excellent reviews by Cowley k Wu 
(1994) and Goldstein (1994) 

Our analysis is similar to that of Wu (1994) who investigated the stability of a two- 
dimensional Stokes-Layer to 3D disturbances and studied the effect of slow- temporal 
and spanwise modulation The major difference m the present work and that of Wu 
(1994) arises from the three-dimensionality of the basic flow used here This leads to 
a different amplitude equation with single spatial derivatives of the amplitude m the 
integro-differential operator and the kernal function, as opposed to the double derivatives 
in Wu’s work. Another important aspect of the current work is the coupling of the 
evolution of the disturbance amplitude with the properties of an unsteady wall layer. This 
extra new feature arises primarily because of the scalings associated with a stationary 7 
cross-flow vortex. 

The basic scaling used in the analysis below may be derived using an argument similar 
to that first used by Hickernell (1984), see also Wu (1994). We piesent some of the details 
of this argument as there are a number of important differences from those given m Wu 
(1994). Consider a cross-flow vortex of amplitude 8. We need to determine the size of 
8 such that critical layer nonlinearity affects the amplitude of the cross-flow vortex. We 
allow for slow temporal and spatial modulation with respective lelative scales A 7 -, A z- 
The suffix Z here denotes variations in the direction normal to the Squire direction Since 
the wavenumbers are of magnitude 0(R * ), R being the Reynolds number, a balance of 
the mertial and viscous terms shows that the thickness of the cutical layer at Y c is given 
by Y — Y c = 0(e) where Y denotes the 0(1) boundary lay 7 er coordinate, and e = R~e. 
Thus the spanwise component of the disturbance velocity 7 , to, has a pole singularity of 
size 

W = °^(Y-Y C )^ (1-1) 

In Wu (1994) the equivalent w is much smaller since 


w 


= 0 ( — £ 2 — ) 


0 ( 


<5 A; 

O'-l'c) 


) 


there. The relation (1 1 ) stems from a different balance Fiom the continuity equation 
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the Squire component of \ elocity u is then of size 


U = O{wz) = Oij0 ^y) 


( 1 - 2 ) 


Again an important difference between the anal} sis here and in Wu (1994) is noted, since 
the equivalent expression theie is 


u - 0 ( — ^ — ) = 0(— Az - ) 
[ (Y-Y C )’ U[ (Y-Y C ) ) 


This is the reason why the final amplitude equation in Wu's work involves double deriva- 
tives in Z whereas here there axe single derivatives m Z 

The remainder of the argument follows closely that given in Hickernell (1984) and Wu 
(1994), and shows that 

A r = A z = 0(Y-Y c ) = 0(t), (1.3) 

for unsteadiness and spanwise modulation to produce a non-equilibrium effect in the 
critical layer. Interactions inside the critical layer give rise to a Squire component of 
velocity of order 

^ £ \ 

l A * T {Y-Y c ) 3h 

and this affects the outer flow provided 


0 { 


S 3 A Z 

A$(Y - Y c ) 3 


) = 0(6A T ). 


Thus using (1.2), (1.3) we find that <5 = 0(e 2 ). 

In section 2 below the problem is formulated The details of the outer inviscid flow 
are considered in section 3 where a solvability condition, which leads to the amplitude 
equation, is derived. The solvability condition shows that the amplitude of the cross-flow 
vortex depends on the displacement induced by the wall layer as well as the jump across 
the critical layer In sections 4 and 5 the solutions inside the critical layer axe obtained 
and these are then used to determine the nonlinear jump conditions which appear in 
the solvability condition. The amplitude equation is obtained in section 5. In section 6 
the properties of this amplitude equation axe discussed and some results are presented 
Finally we conclude with additional comments in section 7. In Appendix A the analysis is 
extended to the obtain the corresponding amplitude equation for a stationary cross-flow 
vortex in the flow over a rotating disk 

Throughout this work the fluid is taken to be incompressible and the Reynolds number 
to be large. 
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2 Problem formulation 


Consider cartesian coordinates (7T, y.z) non-dimensionalised with respect to a lengthscale 
L and where x is m the streamvvise diiection, z is in the spanwise direction and y is normal 
to the body The corresponding non-dimensional velocity components are (u,v,w). It is 
convenient to work in terms of Squire coordinates defined by 

x — ax 4- 0z, z = —f3x 4 az (2 1) 

with corresponding velocities 

u = au + 0w, w = —/3u + aw (2 2) 

where a = cos 0, f3 = sin 6 and the angle 6 will be fixed subsequently. 

In terms of x : z the Navier-Stokes equations are 

div u = 0, (2 3a) 

u ( + u Vu = — Vp + yrV 2 u (2.3b) 

R 

Here u = (u, v,w), t is the non-dimensional time, p is the non-dimensional pressure, 
R = UcL/v is the Reynolds number, Uc is a characteristic velocity scale, and v is 
the kinematic viscosity The Reynolds number R is assumed to be large, and we set 
£ = R~ lf6 . 

2.1 Basic Flow 

The basic flow is taken to be a fully three-dimensional boundary layer flow given by 

{u,v,w) = {Ub,£ 3 V b ,W b ,){x,Y,z) + (2.4a) 

p = p B {x,z) + (2.4b) 

where Y = e~ z y is the boundary layer coordinate. The basic flow satisfies the boundary 
layer equations 

U B x + Vby + W Bz = 0 , 

UbUbx 4 - V b Uby + WbUbz = — PBx + Ubyy- 
U b Wbx + VbWby + W B W B z = ~pbz + Wbyy, 
with the boundary conditions 

U B = V B = W B = 0 on r = o. 

and 

U B - Uboo, W b - W Boo as r — oo. 

Following the discussion m the Introduction, the basic flow U B is such that it has a zero 
at a point of inflexion Thus near Y = Y c where U B {x, 1 c, - ) = 0 we can expand the basic - 
flow quantities as 


U B — eXip , 

(2.5a) 

V B = Vbo + eVgfa + — 

(2.5b) 

W B = Ho + CHiV + £ 2 H2V 2 4 

(2 5c) 
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where 1' = Y c + £y and Aj, A 3 . 1^, \'bo, m 0 ,Mi-/' 2,. and all functions of x and r. 

After substituting (2 5) into the basic flow equations and equating powers of f the 
following relations are easily derived. 


Vbo^I — PoYcz^l — —pBr 

(2.6a) 

V&h, + /^O Air — A \Ycx — Ml Alice = 6A3, 

(2.6b) 

Vbo + MOr — Ycx Al — YczH 1 = 0, 

(2 6c) 

VboPi + MoMor - PoYczPi = ~PBz + 2/i2 

(2 6d) 


These relations are used later in section 4. The properties of the basic flow near the wall 
are also needed later and for Y small we have 


Ub = AiV + (2.7a) 

W B = fhY + . .. (2.7b) 


Since the velocity profile U B is zero at the wall and has a zero somewhere in the flow, it 
is clear that Xi and Ai have opposite signs. 

Next we consider a stationary cross-flow vortex with a wavenumber 7 (scaled with 
respect to boundary-layer thickness), in the x direction, and introduce additional coor- 
dinates X, Z , T to allow for modulation such that, 


d_ 

dx 


£ 


-3 



d 

dx 



(2.8a) 


d _ 2 d d d _ 2 d 
dz £ dZ dz ’ dt dT 

The e -3 factors above account for the scaling ■with respect to boundary layer thickness 
and the J7 variations are needed to account for the local variations of the mean flow 
in the x and z coordinate directions. 



3 Solution outside the critical layers 

3.1 Main part of the boundary layer 

We first consider the solution m the main part of the boundary layer where Y = e~ 3 y is 
0(1) In this region the expansions for the flow quantities are 


u = U B + 

8(u 0 + £U\ + 

) 

(3.1a) 

v = 

6[v 0 + evi + 

) 

(3 lb) 

w = W B -\- 

8(w 0 + ewi + 

•) 

(3 lc) 

P = PB + 

KPo + ePi + 

) 

(3 Id) 


where the disturbance size 8 will be taken to be 0(e 5 ^ 2 ) subsequently, but since the prop- 
erties outside the critical layers are largely dictated by lineal dynamics, it is convenient 
to work with 8 
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After substituting (3 1 ) into the Xavier-Stokes equations 
leading ordei disturbance equations are found to be 

(2 3) and using (2.S), the 

Uoi + vqy — 0 , 

(3 2a) 

UbU 0 £ + v 0 Ub\ = Po£ ’ 

(3 2b) 

Ubw oz + vq Wb) = 0 

(3.2c) 

The problem (3.2) is just that for the stationary cioss-flow 
Walker (1955), and if we set 

vortex, Giegory, Stuait 

= AU 0 (Y, x,z)e ni + c c.., 

(3 3a) 

I>o = —iAVo(Y,x,z)e 11 ^- l-cc.., 

(3.3b) 

= AWo(Y, x, z)e n( + c.c.., 

(3.3c) 

pS = AP 0 (Y,x,z)e n( + c.c., 

(3.3d) 

where c.c.. denotes the complex conjugate, then Vo satisfies Rayleigh’s equation 

'W-(7 2 +% 21 )f'o = 0 , 
Ub 

(3.4a) 

with the boundary conditions 


Vo = 0 on Y = 0, 

(3 4b) 

|V)| <00 as Y — * 00 

(3.4c) 


In (3.3) A(X, T, Z) is the normalised amplitude of the cioss-flow vortex and dependent 
on the slow scales X,T and Z. The eigenfunctions 1 in (3 3) depend on x and z because 
of the basic flow Ub The problem (3.4) has to solved with the additional constraints 


Ub — 0 = Ubyy at y = Ic- (3.4d) 

so that the point Y = Y c is a regular point of (3 4a). The conditions (3.4d) fix the angle 
9 and Y c since 

Ub = cos 9U b 4- sin 9Wb- 

Given 6 and Y c (3 4a,b,c) then determines the wavenumber 7 of the stationary cross-flow 
vortex and the eigenfunction V 0 . and Uq, Wq are obtained fiom 


Uo = -Vo y, 

7 

(3 5a) 

,v° 

II 

Kf 

H— 

(3.5b) 


In general the problem for Vo has to be solved numerically but the solution properties 
near Y = Y c can be obtained using a Frobemus expansion. Thus near Y = Y c the solution 
for Vo can be expressed as 

Vo = C ± <? a + B ± [<j> b + Pj <j> a In \Tj\] (3 6 a) 

x The amplitude A also depends on x and r as can be seen from the amplitude equation in section 5 
below, but since the dependence is only parametric it is not shown explicitly for the sake of clarity 
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where 


9 a = »/ + ~P : V 2 + • , 

<f>b = 1 + + • 


are two linearly independent solutions of Rayleigh's equation and rj = Y — 
Note that 


2Xi 


1 


ba 3 


P) = _=0, = 

The constants 2 C + — C~,B + — B - will later be shown to be equal to zero 
From (3.5b) it follows that Wo has a pole singulanty near Y c if py 0. 
The second order problem for the flow quantities leads to the equations 


+ «o.Y + Viy + Wqz = 0. 


(3.6b) 

(3.6c) 


(3.7a) 


Uox + Ub(uiz + UOA ') + ViUby + Ubuoz = ~P\i - Pox, 

Vqt + Ub{v\Z + Vox) + WbVqz = —Pi) . 

wot + U B {w u + wqx ) + v\Wby + Wbwqz = ~Poz- 


(3.7b) 

(3.7c) 

(3.7d) 


Thus using (3.7) it is found that the e n ^ component of iq. Vj satisfies the equation 


Ub{Viyy — 7 K) — ViUbyy = ~ 


1 t byy -r- dA 


7 U 


B 


T / 

‘ °dT 


1 dA ( Wbyy WbUbyy 


r)A 

UbVo-'I-' — VoUb. 
OX 


7 dZ \ U B Ul 

The solution for Vy near Y = Y c is easily obtained and takes the form 


(3.8) 


Vy ~ 9±<i>a + f± <j>b + 


B ± 1 6 poX 3 2 // 2 ^ dA 6 A3 B* dA 


[ 7 V Af A, J dZ 7 A? dT J 


(77(ln|7?|-l) + ... (3 9) 


Note that since satisfies the same equation as Vo but with a different right hand side, 
a solution exists only provided a certain solvability condition is satisfied The solvability 
condition is obtained by multiplying (3.8) by Vo and integrating from Y = 0 to Y = 00 
and using (3 4) This gives 


where 


[VorK - VoV lY )%± + [WorFiW = + ~ + 


Iy dA I 2 dA 


dA 


7 dT q dZ 


'dX 


h = i 


rUByY VSdY , 


Jo ul 


f 00 (W B yy W b Uby\ \ ,>2 

h = “Jo w~r° dl ' 


(3.10) 

(3 11a) 
(3 lib) 


r 00 _ 

h = -2 / V 0 2 dY (3.11c) 

Jo 

The equation (3 10) is the amplitude equation for the cross-flow vortex and the first 
term involving the jumps is obtained from an analysis of the critical layer The Vy |y =0 
contribution to the amplitude equation is determined from an analysis of the wall layer, 
and this is considered next 

2 These are again functions of x and ; but since the dependence on x ; is parametric we will use the 
term constants 
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3.2 Wall layer analysis 

From (2.7), (3.1), (3 2), it follows that m the wall layei wheie y = s 4 Y the expansions 
are 


u — fAil + . 4- 8(uq + ). (3 12a) 

v = <5(£F 0 + .), (3.12b) 

P = PB + 8(ep 0 + - •)> (3 12c) 

w = efaY + . . + 8(w 0 + - ). (3 12d) 

Substitution into the Naviei-Stokes equations then leads to 

+ v 0 y = 0, (3 13a) 

uqt + Ail uo( + fo-^i = ~ pot 4- uiyy, 0 = —PoY-> (3 13b) 

wot + Ai}'uJ 0 { + Fo/Ji = tF 0 rr- (3.13c) 

These equations have to be solved subject to the no slip conditions 

u = v = w = 0 on T = 0 (3 14a) 


and 

u 0 -» — Vo v(Y = 0)e ; ^ + c c. as Y oo, (3 14b) 

w 0 — » — - bl^l. ZX -ji y QY [Y = 0)e*^ + c c. as Y — * oo. (3.14c) 

7m 

The displacement from the wall layer provides a contnbution to the amplitude equation 
and this is given by the finite part of v 0 as Y — * 00 . Hence the required matching 
condition is 

Vi (Y = 0) = Jim [v 0 + V( ? A(y, Z, T)V qy (Y = 0)e n€ + c.c.)] . (3.14d) 

The solution of the unsteady wall layer equations (3.13) togethei with the non-slip bound- 
ary conditions (3.14a) and slip velocity (3.14b) determine (3.14d). However these equa- 
tions cannot be solved m isolation because of the unknown slip velocity (3.14b) and in fact 
the wall layer equations are directly coupled to the nonlinear evolution of the cross-flow 
vortex amplitude via (3.10). 

We consider next the details of the critical layer which deteimine the unknown jumps 
in (3.10) 


4 Solutions inside the critical layers 

In the critical layer we set y = e 3 l’ c + £ 4 p and the outei solutions, with 8 = £ 5 / 2 , imply 
the expansions 

u = eAjrj -j- c 5/,2 uo + £ 3 (A3^ 3 + 1 / 1 ) + s , ' ,2 x /2 + • • 1 (4 1a) 

V = £ S/2 v_i + £ 3 v$ + e‘ /2 v 0 + e 8/2 (vi + I BO 7 /) + - 9/2£ '2 + •, (4.1b) 

w = Ho + epiy + £ 3/2 u >0 + £ 2 {P 2 T 1 2 + u>i) + t 5/2 it> 2 + . , (4 1c) 

p = Pb + £ 5/2 Po + £' /2 Pi + 1 8/2 P2 + £ 9/2 P3 + • • (4 Id) 
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After substituting (4 1) into the Naviei-Stokes equations and using the relations (2.6) we 
obtain the following sequence of equations: 


Aju_i — — po^, 

(4 2a) 

£■ 

(§1 

II 

o 

(4.2b) 

UO( + Vo v + Wo z = 0, 

(4.3a) 

fCu 0 + r'oAi = — (pi^ 4- Pox)i 

(4 3b) 

ICw 0 + V-ifJ.1 = 0, 

(4.3c) 

0 

II 

»— * 
-a 

(4.3d) 

where the operator K, is defined by 


r _ d . 3 5 52 

L ~dT + lV d£ +flo dZ dv 2 

(4.4) 

The equations for the second and third order disturbance quantities 

are 

ui* 4- v lv 4- wiz = 0, 

(4.5a) 

fCui + tii Ai -f v^iUo v = — P2£- 

(4.5b) 

fCwi + ti-iuith, = 0, 

(4.5c) 

0 
11 

1 

3 

■a 

(4.5d) 

and 


tf2£ + UqX + i’2 JJ + W 2 Z — 1 CzW Ot, = 0, 

(4.6a) 

tCv.2 + klT)U0X + U 2 A 1 + tl-lUlrj + Plljllo z 
+ Ag^liOr, + V-1ZX3T] 2 = — (p3{ + Pi A' )• 

(4.6b) 

Kw 2 + X\T}Wox + U-l Wlv + ^oPi 4- P1VW0Z, 


4- u_ 1 2p 2 r? + V^wor, - PoYczWov = ~Poz, 

(4.6c) 

A.tl_l — P3rj- 

(4.6d) 

The solutions to (4 2)-(4.6) together with matching with the outer flow determine the 
jumps required for the amplitude equation. These equations are solved m senatum using 
the well established Fourier- Transform technique of Hickernell (1984) It is convenient to 
introduce a change of variables with 

Z = p Q 1 Z , T = poT — Z 

(4 7) 

so that the operator K. becomes 


„ _ d . d d 2 
~ dZ + 1?? d(, dv 2 

(4 S) 
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First from (4 2) we obtain 


r_i = -iA{X,T,Z)e ni g x + c c (4 9) 

where we have defined 

A( X, T, Z) = A( X, Ho'T + Z, g. 0 Z) = A{ A\ T. Z), (4 10) 

and gi is a constant. Matching with the outer solution (3 3b), (3 6a) shows that 

9i = B + = B~ (4.11) 

If we put VLi = — tA(X, T, Z)gi, w 0 = W 0 e'^ + c c , u 0 = U 0 e'^ + c.c., v 0 = 
Voe 1 "^ + c.c., then from (4 3) using (4 10), we find that 


W t 


where 


Also 


o = dov f°° A(X, T,Z — s)e'^- ,vni dTi, 
Jo 


1 V HJ 1/Zi 

V - Ai 7, or = — , do = — 

3 Ax7 


U 0 = -—Woz +U 00 (f,Z ), 
ry 


(4.12) 

(4.13) 
(4.14a) 


Vo = -iyn(j 00 + Voo{T, Z), (4.14b) 

where f7 0 o, Voo may be determined through matching In fact from (3.1), (3.3), (3.6), 
(4.14b) we see that 

AC + = AC~ = 7 t/oo (4-15) 

We next consider the solutions of the second ordei problem (4.5). The forcing terms 
m the equation for Wi indicate a solution of the form 

Wi = W 10 + W 12 e 2n * + c c 

where Wio, W 12 are independent of £. From (4 5c) the equation for Wio is 

idw 0 


CoWio = -wii ): 


dg ’ 


where the operator C n is defined by 


„ d n d 2 

£ " = sI + 2m7, -3? 


The equation for Wio may be solved to give 


Wi 


10 


= - g^dov 2 r ds r ^ c \x,T ,Z - Si)A{X,T ,Z - s - s l )sKf ) {s,s l \a)e-"' lls d, 
Jo Jo 

where the notation A^ denotes the complex conjugate of A and 

Af\s u s 2 \a) = e -^ +3s " i2) 


(4.16) 

(4.17) 
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The equation for Hq 2 is 


£ 2 ir 12 = -vi 1 ir ( 


On 


This is solved to obtain 


W l2 

where 

and 


= d x r d Sl [°° siA{ X. T,Z- si)A{X, T,Z- s 1 - s 7 )K l 4 1 \s l . s 2 ) e -*‘ / ’ 7 ( 2si+S2 > ds 2 , 
Jo Jo 

(4 IS) 


K^isuSi) = e-*(*?+ 3i ^+ 6s 2+^). 


d\ = 


(4.19) 


Next if we differentiate (4.5b) with respect to tj and use the continuity equation (4.5a) 
this gives 

ICfli = XlWiZ — V-iflor,, 


and we have defined Cl n = The right hand side of (4.20) suggests that 


fli — fi, 0 -t- fli 2 e 2n ^ + c.c 

The mean flow component of fli, ffxo satisfies the equation 


Coho = + W-4 - 4)11 


10 


dr) ' * vr u dZ dT ' 

Using (4 14a). (4.18) it is found that Clio is given by 

fiio = -g[ c) — V 3 r ds r s 2 a!~ c) (X, T,Z- s 2 )K[ 0 ] (s, s 2 ) X 
7 Jo Jo 

e-'^^o 1 ~ f ' Z ~ " 2 " $) ds 2 

—g^XidoV 2 f ds f ds 2 f sl\[ 0) {s,si + s 2 \cr)e~ un,s x 
Jo Jo Jo 

(fo‘4 - JL)0‘\x,T, Z — So — Sl )A(X,f. Z — So — so - *)] ds , 
The e 2ry ^ component of fli satisfies the equation 

C 2 Cli 2 = — Vliflorj + AiT'Ui 2 ^. 

Thus using the solutions for flo, U / 12 , it is found that 

Ch 2 = -2d 2 v* r ds /’ 0 °s 2 /vi 1) (s.s 1 )e- ,,? " (2si+s) x 
Jo Jo 

A(X,T,Z- S ,)JL(A(X,T,Z - si - s))ds 1 

r oo r oo roo . . 

+4d 3 i/ 4 / ds dsi / Si/v^ (si,s + •s 2 )e _,,,1,(si+2s+2s2 *x 
Jo Jo Jo 


(4.20) 


(4.21) 


(4 22) 
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(4 23) 


d_^ 

dZ 


[A(A\ T. Z — s — s 2 )A(A\ T. Z — $2 ~ *i — *)] ds 2 , 


where the constants are defined bv 


do — 


9\d 0 
22/7 ’ 


do — 


dpffi 

4 2/7 


We turn our attention next to the third order problem defined by (4 6a-d) which deter- 
mines the nonlinear jump. It is only necessary to calculate the e n ^ components of W 2 
and U 2 , W 21 and U 21 respectively From (4.6c) is is seen that 



C\W 2 \ — i?i + R 2 + A3, 

(4.24) 

where 

i?i = — 7?[AiI4oa' + 2^2 F—i + Mi^bc ~ 27/21^00], 
F? 2 = — poz — M1F00 — hfi 0 * TFor; + Ho^czWon-, 

(4.25a) 

(4.25b) 

and 

r 3 = -[\>iWi2, + V-i(Wio, + w 0 (c) »?)] 

(4.25c) 

Writing 

W21 = wff + wg> + h4 3) 


with £iW 2 i 

= R } , then it can be shown that 



w" = - 4^.1 + 

+ - -H'^] + - 4f 0 „„], 

/Ij 22/ 22/ 


where 


-Fo = rUoo(X,T,Z-s)< 
Jo 

Similarly the solution to C\ Wff = R 2 is given by 


—os—ii/ris 


ds. 


MV = — — Woz - -F x + - /2 0 Vcc)ir 0w , 


where 


(4.26) 


HMl Mo 

Fi = r Voo(X, T,Z — s)e~° s3 - iur,s ds 
Jo 

Next, using (4 9), (4.16), (4 18), the solution to C\W^ = i? 3 is found to be 

wg ] = -d 4 2/ 3 r ds 2 r ds 1 [°° s 2 2 A{X, T,Z-s + s 2 )A {c \x. f.Z-s + s 2 - Sl ): 
Jo Jo J $2 

A{X, T,Z — s — si)e~ u/STI Sl .s 2 ) ds 

+d 5 1 / 3 / ds 2 / dsj / s\A(X,T. Z — s — s 2 )y. 

JO JO J — 52 


(4.27) 


A(X,T,Z- s-s 2 - s 1 )A [c) {X,T,Z - s - 2 a 2 - s 2 ) ds 


— ifs 7 ) r,'( 3 )/ 
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A{X, t , Z 
where 


-del' 4 r ds 2 r ds r s 25l ,4 (c) (A'. T Z - s + < 2 ): 
0 J S 2 J 0 


1(25 - 5 2 - Si))i4(-Y. T,Z- ^(2s + Si - s 2 )K^ ) {s.s 1 .s 2 )e- ,,n >‘ ds u (4 28) 


Ki 2 \s,s u s 2 ) = e-^ S ^K 

/ii 3) (5,5!,S 2 ) = e-A^l+Sshi+s 3 ) 


d 4 = |<7i| 2 do, ds = g\d\^\ de = . 

iv 


(4.29a) 

(4.29b) 

(4.29c) 


Finally, one other quantity which is required for the calculation of the jump is f I 2 i, the 
e n * component of fi 2 which satisfies 


C 1&21 — Qi + Qi 4- Q 3 + Qii 


with 

Qi = — r](XiClox + 6 A 3 VI 1 + fiiCloz + 7 2 A!F_i). 

Q 2 — —pii-oz + tyV_iz — Vbo&ov + AilcsH'oij. 

Q 3 = -V z i(Cli Or, + ft&M - 
Q 4 = \iW 2 iz- 


If we write fl 2 1 = fI 2 V + fl 2 i^ + & 2 i + f^ 2 ^ with then it is eas}' to show that 


fi&> = ^ 
V mi 


Ifl 


/Xl a iivWom-^Wonnm). 


„ d 2 

+ — ( 

dZ 2 

The solution to = Q 2 is using (4.3c), (4.12), (4.14). (4.26), 

dF 0 ndW 0 , tf, Ail )-- 7 ~ 7 

'~dz ~ 7 i ~az + ~ -*T 1W 

Next &21 satisfies 


where 


r\(2) Ml d 2 , l3 , , 3F 0 ?7 dWo , VJ 

= 7a7 ( ~ H lz~ 77z + 7 

ZidS = q? + q? + q ^ 3) - 


Qi x) = -V-A 0 ,, Qi 2) = -V-I«a, Qi 3) = -ri?Qi 2 , 


The solution Cl^ is decomposed into 


ng = E 

j=i 


(4.30) 


(4.31) 
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with -Cifi^i'"^ — Qz \ an d each of the O 2 I J ^ components aie sohed foi separately However 
only the ^ teim piovides a non-zeio contubution to the amplitude equation Using 
(4.22), it is found that 

^ 2 i’ 2) = -del' 4 r ds * r ds * r 4A(X. T,Z-s- s 3 )x 
JO Jo J — S 3 

A(X, f , Z — s 2 — s — s 3 )I\[ 3 \s, s 2 , s 3 )e 11/7,5 ,f. Z — s 2 — s — 2s 3 )j ds 

+dju / ds 3 ds 2 / ds\ / slA(X,f,Z-s — s z )I\l\s,s 2 + Si,s 3 )x 

Jo Jo Jo J—sz 

e-'^-^lAiX, f,Z-s-s 3 -s 2 - 5l )I (c) (X, f , Z - s 2 - ^ - s - 2s 3 )] ds, (4 32) 


where 


de = _sML, dr=ti4” V 
7 


Finally the equation for Oji 1S 

= a 1 ww = £a 1 iv 2 ( & 

Writing ft 2 ^ = E;=i ^ 21 ’^ where 


3 

z 

J=1 


Ang 1 ’’ = a,h' 2 “ 


(4.33) 


from (4.26), (4.29a) we obtain 


^ _ + 2 7 2 aza * [r)Wom 2iv 617V azax h/oi?W! 


1 a 2 


1 


1 a 2 


mi a 


1 


a 2 


Ml t/ r , _ _Lu/ i Ml t/ s 

2 7 2 A 1 aZ 2 ^ ^ 2iv 0rmm] 6i lt AdzJ 0vvnr,) 

~ + b*z w °~ 1 


Mi 


a 2 1 a 

— ^TT7^ 71 ° ,7 ’ 7T7 -l + 


3 zj/ 


3 z/i az 


(4.34) 


and 


^ + ^Jz F " + ^ (vi » - * 1 ->S2 (,,r — >■ (4 - 35) 

In (4.28) the solution for imolves three triple integrals and only the second triple 
integral generates a term m the solution for S2 2 l’ 3 ^ which contributes to the amplitude 
equation The solution to (4 33) with just this term present is 

^ 2 i’ 3) = Aid 5 z/ 3 f ds 2 [ ds f (s 2 + s)sll\l t 3) (s.Si,S 2 )e~“ /s ’ , x 
Jo Jo J —$2 

^r[A(X, T,Z — s — s 2 )A(X,T , Z - s-s 2 - s^ c \ X, T.Z-s- 2s 2 - s 2 )] ds. (4.36) 
We turn our attention next to the calculation of jumps 111 (3.10) 
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5 Calculation of the jumps and the amplitude equa- 
tion 


From (3.1), (3 9), (4 lb), (4.11), (4.14), (4.16) it is cleai that (\\)\ e c t = 0 In addition 
these equations show that 


(Kr)k + = (g+-g-) = -t 7 / (fi 2 i )cty, 

Jt )=— oo 


(5.1) 


where the notation /* denotes the finite part of the integral The integral m (5.1) is now 
evaluated using the solution for 0 2 i obtained in the pievious section. First from (4.30) 
we find that 

H = A (6 a 3 + 7 *A 1 K«43(A- 1 T, Z), (5.2a) 

j — oo oZi 

where s„ = sgn(Ai), and only the first term in (4.30) gives a nonzero contribution. 

Next using (4 31) we have 


r^dr, = -s n ^^(X,f,Z)-s n ^^-^A(X,T,Z). 

J - oo v oZ /ii dZ 


d 0 7T7 7 8 , - - , 


(5.2b) 


Whereas (5.2a, b) are linear contributions, the f^i’ 2 ^ terms gives a nonlinear jump From 
(4.32) we find that 


/ *oo _ . x roo roo 

Clfidrj = —s n deu 3 2-ir ds 3 slA{X,T,Z — s 3 ) 

•oo */0 Jo 


X 


A(X, T,Z — s 2 — s 3 )hf\0, s 2 ,s 3 )^ c> ( A', T.Z- s 2 - 2 s 3 )] ds 2 

roo roo roo _ 

+2s n 7rd 7 i/ 2 / ds 3 / ds 2 / ^^(A, T\ Z - s 3 ).I\ 4 '(0, s 2 + si, s 3 ) 

jo Jo Jo 

*-\ 

^[A(X, T , Z — s 3 — s 2 — Sl 0 c) (X, T,Z- s 2 - Sl - 2s 3 )] <fa,. 
oZ 


& rrM/ 


Next the contribution from the term i 


is 


2 M2 do* d t/ v « 


Vi* d 


r ^' 1)dT l = Sn m2=-j=A(X t T, Z) + 5»— ^^0 °(A', f , Z). 
j-oo z-y/zi c/Z */ aZ 


(5.2c) 


(5.2d) 


The f^i’ 2 * term gives zero contribution. The f^i’ 3 * term gives a nonlinear jump term 
f ni 4 { 3) dr) = s n A 1 d 5 i' 2 2r [ ds 2 [ dsi sfA'k^O, si, s 2 ) x 

J— oo Jo Jo 

^[A(X, T, Z - 5 2 )A(X, T , Z — s 2 — Sl )A {c) (X, T.Z- ls 2 - Sl )] (5.2e) 

Hence finally collecting all the terms (5 2a-e), restoring the variables X, Z, T, and sub- 
stituting into the (3.10) gives the resulting amplitude equation as 
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h dA I 2 dA dA 
jdT + -fdZ + lh dX 




f 9\* /0 6A 3 /x 0 &4 gi a6A 3 dA 

AiSnigi i — -(2fi2 ; — - 


A l7 2V "^ Ai 'dZ A 2 - 2 dTj 

K CO TOO 

*3 ^2 5 3 A(A'’, T — ^3, Z — //oS 3 )A ' 4 (0,62153) x 

A(A, T - S 2 - S 3 , Z- Ho(s 2 + 5 3 ))^-^-/i^ C ^(A, T — S 2 — 2 s 3 , Z — fJ.o(s 2 + 2 s 3 )) 

roc roo roc 

+ ds 3 / ds 2 / dsx 5|A(A,T- 6 3 ,Z - /Zo53)A'i ) (0,5 2 + 5l,53)x 

Jo Jo Jo 

^:[A(X, r-S 3 -S2-Sl,Z-//o(63 + 52 + Si ) A^ (A, T — S 2 ~ 61 — 253 , Z — /Zq(S2 + 61 + 2S3)] 


+ 


r ds 2 r 

Jo Jo 


5 

ds! SjA'f^O, Si ,s 2 )-^[A(X,T - s 2 . Z - fi 0 s 2 )x 


A(X , T — S 2 — S\, Z — flo(s 2 + Si))A^(A, T — 252 — Si, Z — Ho(2s2 + 5i)) ] (5 3) 

This can be written in a more compact form as 


hdA I 2 dA dA - 

t"5T + ~^ydZ ^^ 3 dX ~ + M(Ji + J 2 + J 3 ) 


(5.4) 


where 


M = 2s n ir'r 2 gl\gl\fi 1 Xl, 



Uby 


Yr~j 


Ul 


V 2 dY, I 2 = - 


B 



WbUbyy \ 
U 2 b ) 


V 2 dY, 


roo roo 

Ji= ds 3 / ds 2 slA(X, T - s z ,Z - iio$z)Kf ] {s 2 , 5 3 |a) x 

jo Jo 

o 

A( X, T — 52 — S 3 , Z — go(s 2 + S 3 ))-^A^(X, T — s 2 — 2 s 3 . Z — go{s 2 + 25 3 )), 

roo roo 

J 2 = / ds 3 / ds 2 s\s 2 A(X , T - s 3 ,Z - ^s 3 )I\^\s 2 , s z \a) x 

Jo Jo 

o 

— — [A(A, T — S 3 — 52 , Z — go {$3 + S 2 ))A^(X, T — s 2 — 25 3 . Z — go(s 2 + 253 )], 
oZ 

too roo . , Q 

J3 ~ Jo dS2 Jo dSl S * K ( S i’' S2 l Cr )aZ ^( A - r “ ^ 2 -Z - g 0 s 2 )x 
A(A, T — s 2 — Sj, Z — go(s 2 + s i))A^(A, T — 2s2 — -si, Z — /zo(2s2 + Si)) ] , (5.5) 

and the path of integration C for A, I 2 is along the positive real axis with an indentation 
below/above the point Y = Y c depending on whether Aj is positive/negative. 

Equation (5 4) is the mam result of this paper. A generalisation of this equation for 
the flow over a rotating disk is given m Appendix A. In the next section we consider a 
few of the properties of this equation and the numerical solution of a special case. 
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6 Special cases and results 


The solution of the full equation (5.4) is m general quite difficult because of the cou- 
pling with the wall layer and secondly because of the form of the nonlinearity. There 
are however some special cases which can be considered further Below we have used the 
normalisation Voy (0) = 1 for the solution of the leading order eigenvalue problem (3 4 ). 


(a) If we consider the plane wave A = Ae^ aiX+/3,z ^~ ,u,lT with Qi./?i,u 7 all real, then the 
contribution from the nonlinear terms J\ + J 2 + Jz is identically zero and (5.4) leads to 
the linear dispersion relation for neutral waves as 


7 


l hP 1 T 

1- 17/30:1 


lt X 6 *" 

(ri^p 


G(£o), 


( 6 . 1 ) 


where 


£(£ 0 ) — £0 + 




{0 = 


(tI^i 1)3 


The right hand side of (6.1) is obtained from the solution of the wall layer problem (3.13). 
The real and imaginary parts of (6.1) gives two equations which can be solved to obtain 
correction terms to the wavenumber and wave-angle for a given frequency 07 . With 07 set 
to zero in (6.1) we obtain the linear neutral results of Hall (19S5) for stationary cross-flow 
vortices. The case u> 1 ^ 0 case therefore is a generalisation of the Hall (1985) results to 
‘almost’ stationary vortices. The function G(£ 0 ) in ( 6 . 1 ) was calculated numerically and 
is shown in Figure 1. For large |u;i| it is easy to show that 


G(( 0 ) ~ 

G(£ 0 ) ~ 


(oe ‘« 5n (7lAi])3 

51^ n <fo(7lAi|)3e 1Sn « 
72z Ja-*! 1 2 


$„uJ 1 — + OO 


as SnCJi — > — 00 


This implies that the influence of the wall layer on the cioss-flow vortex diminishes as 
the scaled frequency increases. 


(b) If we consider a disturbance of the form A = A(T)e’^ aiX+/3lZ ,wlT ) with real 

and ai complex, substitution into (6 1 ) shows that 


— ta 1 


Ujli fi\l2 

[-1 \-i 

7 7 


+ i 


e'S’-Gito) 1 

(l|A,|)i ih 


(6 2) 


Since / 1 , 12 are complex it can be seen that disturbances of this type are unstable with 
the growth rate increasing indefinitely for -large |/3i| 01 |_'i | 


(c) Consider a disturbance of the form A = A(Z)e'^ aiX ~ UJlT \ where 07,07 are real. This 
form of disturbance could, for instance, represent a combination of waves of the type 
considered in section (a) above. In addition this is a special case relevant to rotating disk 


IS 



flow, (foi example if we take B = B{R)t'^ @ u ' lT) , with leal, in (A6) of Appendix 

A) The equation (6 1) then leduces to 




(6 3) 


where 


J = I jf ds 3 J^ ds 2 slK ( 4 °\s 2 ,s 3 \cr)A{Z - s 3 )A(Z - s 2 - s 3 )-^-(Z - s 2 -2s 3 ) 

+ / ds 3 [ ds 2 sls 2 I\[ 0) (s 2 , s 3 \d)A{Z - s 3 )AL[A(Z - s 3 - s 2 )A {c] (Z - s 2 - 2s 3 )] 

Jo Jo aZ 

+ 1 ds ^J Q ds 2 sll<i°\s 2 ,s 3 \d)-^[A(Z - s 3 )A{Z - s 2 - s 3 )A {c) (Z - s 2 - 2s 3 )], (6 4) 

and we have set 


k = 


7_ 

h 




7 


— ux\il 3 — i- 


■<?(&)' 


(7|Ai|)i 


i 

K 


h\no\ 5 - _ cr 
M ~ |^o! 3 ' 


(6 5) 


As m Goldstein k Leib (1989) it is found that the numerical solution of (6.3) points to 
a singularity a s Z ^ Z s A local asymptotic description of the singularity can be found 
by writing 


( Zs-zy 


( 6 . 6 ) 


with a,r complex constants When (6 6) is substituted into (6.3) it is seen that for a 
balance of the dominant terms on the left and right hand sides of the equation we require 
t = | + iTo, where tq is real. This leads to 


J ~ 


(z, - zy» 


D{r b) 


(6.7) 


where the function D(t 0 ) is given m Appendix B and D(r 0 ) satisfies the equation 


D(tq) _ (j> 

t |a| 2 /c 


(6.8) 


In Figures 2,3 we show r 0 and |a| as computed using (6 8) and (Bl) as a function of 
- ar s(£) = ar g( C -~ ^ ) In Figures 2,3 we have taken |^>| = [k| = 1 since as in Goldstein 
k Leib (1989) equation (6 3) is completely characterised by the arguments of 6 and k 


6.1 Numerical Solution of equation (6.3) 

The equation (6 3) was solved with the initial condition 

A — * as Z — * — oo. 

Here we have assumed that /j. 0 is positive The form of the arguments m the kernal 
function would suggest an apparent difficulty when Ho is negative It is seen that when 
Ho is negative the nonlinear terms in equation (6 3), (6 4) imply that the amplitude at 
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the current Z location is calculated from a knowledge of ,4 at positions ahead of Z , in 
contrast to the case when /<o is position when the calculation of A(Z) involves history 
dependent effects This apparent difficulty can howevei be lesolved by redefining 9 used 
in (2 2), (by for example adding a multiple of ±tt as this does not affect the subsequent 
criteria used to fix 6) Similar comments apply also to the lotatmg disk flow and the 
appropriate choice for fi 0 there, see (A6), (A7), is to take /i 0 to be negative. 

Equation (6.3) was solved numerically using a 5th oidei Adams Bashforth predictor- 
corrector scheme. The integrals were truncated to those over a finite domain and eval- 
uated using a trapezoidal rule One novel feature of the present implementation of the 
method worth mentioning is that the computations were performed on a massively par- 
allel computer, the Maspar MP-1104 which has 4096 processois arranged in a (64x64) 
square matrix. The bulk of the computational time m solving (6.3) numerically arises 
from the evaluation of the integrals especially when the nonlinear terms become signif- 
icant. Integrals of the type occurring in (6.3) can be evaluated extremely efficiently in 
parallel. Further details of the algorithm used may be obtained from the author. 

Some solutions of (6.3) obtained numerically are shown in Figures 4-6. In Figure 
4(a-c) results are presented for the case with arg(/c) = 0 and arg(©) = — 7t/2 for values 
of a = 0, 0.5, and 5 In these figures the dashed line is obtained from the asymptotic 
solution (6.6). The computations show that the location of the singularity is delayed 
with increasing a. As the singular point is approached there is a sharp reduction in 
the wavenumber, Figure 4(b), accompanied by a very laige increase in the growth rate, 
Figure 4(c). The main solution characteristics for this set of parameters is broadly in 
line with those found by other investigators in their studies of related mtegro-differential 
equations. 

The solutions presented in Figures 5(a-f) , with arg(«-) = 0,arg(<^) = — 7r/4 and a = 
0, 0.5 and 5, however, show a number of new and mteiesting pioperties, some of which 
have not been found before. The results for cr = 0 and 0.5 are similar to those in Figure 
4(a-c) and show again that the singularity is delayed with increasing a. For a = 5 on the 
other hand, our results, up to the largest Z value that we have been able to compute, 
indicate that the singularity has been eliminated in favour of a large amplitude nonlinear 
oscillation. Figure 5d shows that the wavenumber fluctuates and contains a large high 
frequency component causing the wavenumber to peak at specific locations. The growth 
rate is seen in Figures 5(a), 5(e), to oscillate about zero, with again very large peaks 
near specific locations. In Figure 5(f) we show \A\ against Z and this shows cleaxly the 
development of the nonlinear oscillations after an initially exponentially growing linear 
phase. This type of solution has not been found in studies of other related integro- 
differential equations. 

In Figure 6(a-c) we present results for the case arg(/c) = 0, and axg (<f>) = 7r/4 for cr = 
1 and 5. The solution properties are broadly similar to those for the case with <f> = — 7r/4, 
except that the wavenumber has large negative peaks at certain locations. The real part 
of the growth rate for a = 5 is exactly the same as that in Figure 5(e). 

The comparisons between the asymptotic and the numerical results in Figures 4-6 are 
quite good and this indicates that the correct singularity structure has been captured. 
In plotting the asymptotic predictions, for a given value of aig(T>(r 0 )/r) the value of r 0 
and a were obtained fiom (6 3), with the value of Z s extrapolated from the numerical 
results 

The numerical solutions take an extiemely long time to compute, expecially when the 
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nonlinear terms become significant Fuither more extensne calculations are currently m 
progress to exploie a wider range of parameter values, and to see whether the nonlinear 
waveform foim found for some parameter values, persists 01 is damped out for large Z 


7 Further discussion and conclusions 


In this paper we have obtained a novel mtegro-differential equation which describes the 
nonlinear evolution of stationary cross-flow vortices in thiee-dimensional incompressible 
boundary layer flows. It has been shown that the evolution of the vortex depends crucially 
on the dynamics of the unsteady critical layer as well as the dynamics of an unsteady 
wall layer. In companion papers the work presented here is extended to non-stationary 
cross-flow vortices, Gajjar &: Arebi (1985), and to compiessible flows Gajjar &; Sibanda 
(1985). In Gajjar & Arebi (1985) it is found that the amplitude of the non-stationary 
vortex satisfies a similar equation but without the wall coupling present. In addition this 
equation has an additional Hickernell (1984) type term whose coefficient depends on the 
curvature of the effective velocity profile at the critical layer. The influence of the wall 
layer in the current problem diminishes as the the magnitude of the scaled frequency 
increases as was shown in the previous section. 

The amplitude equation has a number of interesting pioperties some of which have 
been discussed already. The full problem (5.4) is of considerable interest and merits 
further study both analytically and numerically. As in many related problems it has been 
shown that solutions to the amplitude equation can develop finite-distance singularities. 
A preliminary analysis of (5.4) suggests that focussing type singularities of the form 


A(X,Z,T) ~e- ,WlT 


F{X) 

(Z s -Z) l + ' T0 ’ 


A' = 


(A- - X.) 

(Zs-Z)' 


may also exist. The function F(X) satisfies a nonlinear first order integro-differential 
equation which can be written down, see also Wu(1994) 

The coupling with the wall layer found here is impoitant m another context, namely 
the study of the receptivity of stationary cross-flow vortices to suiface mounted obstacles. 
Experimentally it has been observed that even minute roughness elements can act as a 
trigger for stationary vortices, see Wilkinson et al (1983), Reed & Saric (1989). It is 
suggested that the close coupling with the wall layer and the manner in which this affects 
the evolution of a cross-flow vortex, may m fact provide a simple explanation for this 
phenomenon. The scales and structure presented here may be used to study this aspect 
in more detail. 

Solutions of the amplitude equation m which the finite-distance, or focussing type, 
singularities form, although mathematically interesting do not however tie in with the 
observations m some experiments and numerical simulations of a nonlinear saturation 
of stationary and non-stationary cross-flow vortices In this repect some of the other 
solutions shown m the previous section, in which a nonlinear wave develops, may have 
more relevance. Other possible equilibrium solutions of the equation are currently being 
investigated. With other scalings, see Gajjar (1994), it has been shown that the evolution 
of long wavelength cross-flow vortices is governed by the full unsteady nonlinear critical 
layers equations In many related problems where similai equations arise, it is typically 
found that the growth rate of the disturbances is driven to zero Thus this type of 
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critical layer nonlinearity may also piovide an explanation foi the nonlinear saturation 
of the vortices 

Stuart in GSW (1955) found that the number of voitices predicted by the linear 
mviscid theory was much greater than that observed m then experiments. Although a 
number of suggestions have been made to account for this disciepancy, our computations 
show that nonlinearity provides a wavelength mcreasing/deci easing mechanism. On the 
other hand the flow in the neighbourhood of the singularity, where the wavenumber 
is changed by an 0(1) amount, is no longer governed by linear dynamics, but rather 
the full Euler equations. A detailed comparison with experimental and other data is 
clearly desirable but requires substantial further work A simple evaluation of some 
of the constants arising from the linear mviscid eigenvalue problem (3.4) is clearly not 
sufficient. 
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Appendix A 

In this appendix we considei the flow ovei a lotating disk and dem e the corresponding 
amplitude equation for a stationary cross-flow vortex 

Consider a disk which rotates about the r axis with angulai velocity Cl. Relative 
to cylindrical polar corrdinates (r, 6. r) which rotate with the disk, the continuity and 
Navier-Stokes equations, suitably non-dimensionahsed. aie 


V u = 0, (Ala) 

^U. 1 

“ + (u V)u + 2(k x u) - rr = -Vp + -V 2 u, (Alb) 

Ol jTL 

where u = ( u,w,v ) are the velocity components, p is the pressure, r and k are unit 
vectors in the i — and z— coordinate directions, and R is the Reynolds number. The 
Reynolds number is taken to be large. With e = i? - s and r = e 3 Y. the basic flow 
is given by Von-Karman’s exact solution of the Navier-Stokes equation, (u,w, v) = 
(u(Y),rw(Y),e 3 v(Y)),p = p{Y) where u,w,v,p satisfy 

u 2 — (1 + w) 2 + uv = u", 2u(l + w) + iu'v = w", 

v + 2u = 0, p' + vv' — v" = 0 (A2) 

The boundary conditions for ( u,w,v ) are 


u = w = v = 0 on Y — 0, 
u = 0, w — » — 1 as Y — *• oo 

Next consider a stationary cross-flow vortex at a location (r 0 , 6) and introduce a 
multiple-scaling as in section 2 such that 


d ,f a d 

dr^ C [ a °d£ + e d£. 



d 

de 


.-3 


' d d 


d_ _ 2 d_ 
dt ^ £ dT ' 


Here ao = cos Oo, 0o/ro = sin Oo, and 0 o is chosen such that (qoU + (&ol r o)w) has a zero 
at a point of inflexion This also fixes the location of the cntical level Y = Y c The 
expansions for the flow quantities are similar to (3.1) with 


u = u + 6(u 0 + eui + . ), 
w = rw + 6(wo + ew\ + . ), 

v = t 3 v + d(u 0 + evi + . ), (A3) 

P = P+ <HPo + e Pi + )- 


After substituting (A3) into (Al) we obtain 


= , Po - = n 

O 0 U 04 H + Uo y — 0 , 

ro 

(A4a) 

Ct 0 u u + uor "1 -j w oe + i’ll' = 0 

^0 r 0 

(A4b) 
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(a 0 f< 4 —U')vqz 4 v 0 Ti) = — a 0 p 0 £- 
To 

, _ , = 

(a 0 u H + u 0 uiy = p 0 £- 

J'o » o 

(a 0 u 4 —w)voi = -Pov, 

To 

Sot + (QoS 4- —w)u 1£ 4 wu 0 h + — u 0 @ 4 i'iUy = -QoP u ~ Pori 
To r o 

- . ( - , A)-»- ,= - /^o -=. 1 = 

to'or + (a 0 tt -I w)wi f + uw 0 R 4 i^oe 4 t’i «■'} = p 1£ p oe , 

r 0 r 0 r 0 r 0 

= /?o _ __ ir _ 

u OT 4 (£*0U 4 4 Ul’oH 4 t’o©4 = ~PlY 

r o t 0 

Next if we define 


(A4c) 

(A4d) 

(A4e) 

(A4f) 

(A4g) 

(A4h) 

(A5) 


-8 0 

X = cxqR 4- ft©, % = -R 4" ct oTo0, 

t 0 

17b = a 0 u 4 — ti>, Wb = — — u 4- Qod). 
t 0 t 0 

= , ft = _ 0 o = 

= a 0 Ufc H ujfc, in* = UA 4 a 0 U'A, 

To To 

then using (A4), (A5) it is found that the equations foi (ii 0 . u 0 .u>o), (fii, tH,t&i) are iden- 
tical to (3.2) and (3.7). We can thus introduce a similar normal mode decomposition so 
that for example, 

u 0 = A{X,T,Z)U 0 (Y)e^ + cc 

The expansions for the wall layer and critical layer follow analogously. Thus in the critical 
layer where z — e 3 Y c 4- eft we have 


u = n 4 e^tio 4 efti 4 e 5 t £2 4 e J U3 4 e 2 u,j 4 - - • , 


3 - 


w = rw 4 e 2 u>o 4 e 2 u>] 4 e 2 u >2 4 4 e 2 4 


5 A 


,3 * 




1 . 


v = e°u 4 e 2 C-i 4 4e 2 i)o 4 e 4 ui 4 e 2 f ’2 4 ... 
p = p 4 e 2 p 0 4 e 3 pi 4 e^P2 + e 4 p 3 4 e^p 4 4 - - 


If we define 

A* = 

and 


d k 


ft d k _ 


-o sfrW + ^31* 


Pfc = 


y=y c 


ft _ </* 

i 0 dr iU + Qo dr fcU; 


V=Ve 


C7t = aotifc 4 —Wk, I V k = ~—Uk 4 a 0 Wk . V k = v k . P k = p k 

TO T 0 

then to the required older ( U k , V*, Ifft -Pfc) satisfy the same equations as the corresponding 
variables in section 4 The amplitude equation, in terms of A. is therefore identical to 
(5.4). In terms of R, 0 variables however, if we define the scaled amplitude of the cross- 
flow vortex as B(R,T.Q ) where 


B{R, r,0) = A(X,T.Z) 
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and use the definitions of A’ Z fiom (A5), then the amplitude equation is given by 
hdB I 2 ( fodB aodB\ ( dB do dB\ .... .... . 

^ /3 V 0 ^ _[ oi lh=0 ( 1 2 3) ' 

(A6) 


where 


Ji = /°°d5 3 r ds 2 slB{R+^^,T-S3,e-^^)Kt\s 2 ,s 3 \(7)x 
Jo Jo r 0 ?'o 

D/'D , MoA), , \ Tn ^ Qo (52 + 5 3 ) ( A> d a 0 d \ 

B(R + — (* + »3 ), T - ,, - ,3, e - «>■——) + -ae J - 

fsW(« + ^(Sj + 2 s 3 ). r - 32 - 2 s 3 , 0 - —(32 + 2*>))1 , 

L r 0 I’O 

h = r ds 3 r ds * 4 S 2 B(f2 + —53, r - S3, 0 - — 5 3 )Af ^2, 5 3 |flr) X 

JO JO To To 


JL i Q ° ^ \ 

v r 0 5.R r 0 50 y [ 


£( A + —(52 + 5 3 ), T- 6 3 - 32, 0 - ^(53 + 52)) 
t~o r 0 


B^( R + ^(s 2 + 2s 3 ), T - a 2 - 25 3 , 0 - f —(s 2 + 2s 3 )) 

1~0 r 0 

* = r *■ r *• w HNs + sNs) [*<* + 


B(R+ ~———(s 2 + 5i),T-s 2 -5 l5 0 - ^^(52 + 5l))x 
7*0 7*0 

B^(R + ^(252 + 5i), T - 252 - 5a, 0 - — — — — (2s 2 + 5a)) 

T*/n 7’a 


(A7) 
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Appendix B 

The function D(r 0 ) is defined as follows with r = j + it 0 


where 


D(t 0 ) — — 2t^I\.i + t(I\ 21 + A22) — ^ ( A 23 + A24) — "(A31 + A32), (Bl) 

a'i=- r ip r«- {p ~ 1] " 

jo jp p 1 q T 

roo roo 

A '21 = / dp dq—— 

Jo Jp P 0 w , 


p p T q T {p + q - l) T<cl + 1 ’ 

03 /-os (p — l) 2 

P p r q T {p + q- 1) tU) ' 

roo rOO 


roo roo 

#23 = / dp dq —— 

Jo Jp p T q T 1 


p T-1 § T+1 (p + 9 — 1 ) t(c) ’ 

(P-1) 2 


L24 


p P T 9 T-1 (P + 9- 1) t(c)+1 ’ 

00 /-oo (p-1) 2 

a( K7- —- 

V 

00 roo 


roo rOO 

“ Jo d? Jp dq p T ~ l q T ( 


p T ~ 1 q T (p + q - 1) T<C)+1 

# 32 = f°° dp I" dq- (P_1)3 
Jo Jp p T 


'p P T 9 T+1 (P + 9 — 1) tU> 

These integrals can be evaluated as in Goldstein &. Leib(19S9) to give 


#1 = E 


(_l)(n+l) 1 


n(n+ 1) (^ + iT 0 ) r 


113 
(-- - ?7 0 ) n - 4(- - ?T 0 ) n + 6 (- ~ 7T 0 )„ 


n «- 


5 7 1 

-4(- - ?To)„ + (- - 7T 0 ) n j , 

#21 = A ( a 2) (r 0 ) - 3A (2 ) (t 0 ), 

2 


#22 — X) 


(_l)(n+D 


^n(n + l)(n-l)(§ + ir 0 )„ L 2 


(-0 “ 7T °)" “ 4 (~A “ 2T °)n + 6(i - lT 0 ) n 


O J 

— 4(- — 7T 0 )„ + (~ — lT 0 )n 

.. ^ (“IT 1 T/l x o,3 . „,5 , ,7 . 1 

23 = 5 »(" + 1) (| + ,r 0 )„ l ( 2 “ ~ 3( 2 " ' To) " + 3( 2 " ’ T ° ) " “ ( 2 ” "°H 


A 24 — X) 


(_l)(n+l) 


^2 n ( n “ l)( n + 1) (f + lT 0 ) n 


113 

(— o ~ 7T o)n — 4(- — lTo)n+ 6(- — lT 0 ) n 


5 7 

— ^(9 — **<,)» + (— — ?T 0 )„ 


# 31 = Ai 2 ) (r 0 ) - 3Ai 2) (ro), 
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a' 3 2 - Y2 


(-i) n 


^ n(n + l)(j + ir 0 ) n 


(-9 “ tT 0)n ~ 4(-~ - l7o)„+ 6 (- - lT 0 ) n 


-4(“ - lTo) n + (- - ITq) t 


where 


1X1 1 

A ‘ I,(T0) = 5 n(n + 1) (| + ~ito\ 1(1 - r - ’ ro) '‘ - (1 + r - ir ° ) ” 1 ’ 

Apl(ro) = S B ( (n+i)(i+.n,)(i+^j : 1(1 - r - !7o) " - (1 + r - * To)J ’ 

and (a) n denotes the function F(a + n)/T(a). 
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Figure 1. — (a) Real (solid line) and imaginary (dashed line) 
parts of G{e~ l against p. (b) Real (solid line) and 
imaginary (dashed line) parts of G(e -# ^ 6p ) against p. 
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Z 




Z 

Figure 4.— (a) A plot of log |A| against Z as calulated numerically (solid line) 
from the solution of (6.3) and from the asymptotic solution (dashed line) with 
arg(K) = 0, arg(<J>) = -it/ 2 The labels (i) and (ii) on the graphs are for 5 = 0, and 
5 respectively, (b) Imaginary part of A'/A against Z, other parameters are as in 
Figure 4(a). (c) Real part of A'/A against Z, other parameters are as in Figure 4(a). 
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z 


Figure 5. — (a) A plot of log |A| against Z as calculated numencally (solid line) from the solution of (6.3) and 
from the asymptotic solution (dashed line) with arg(K) = 0 arg(4>) =_-ti/4. The labels (i), (li) and (iii) on the 
graphs are for 5 = 0, 0.5 and 5 respectively, (b) Imaginary part of A'/A aganist Z, other parameters are as 
in Figures 5(a). (c) Real part of A'/A against Z, other parameters are in Figure 5(a). (d) Imaginary part of A'/A 
calculated numencally with 5 = 5, other parameters are as in Figure 5(a). (e) Real part of A'/A calculated 
numencally with 5 = 5, other parameters are as in Figure 5(a). (f) A plot of |A| as calculated numencally 
with 5 = 5, other parameters are as in Figure 5(a). 
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Figure 5. — Continued, (c) Real part aiA'/A against Z, other_ 
parameters are as in Figure 5(a). (d) Imaginary part ofA'/A 
calulated numerically with 5 = 5, other parameters are 
as in Figure 5(a). 
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Figure 5. — Concluded, (e) Real part ofA'/A calculated 
numerically with 5 = 5, other parameters are as in 
Figure 5(a). (f) A plot of |A| as calculated numencally 
with 5 = 5, other parameters are as in Figure 5(a) 
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Figure 6.— (a) A plot of log |A| against Z as calculated 
numencally (solid line) from the solution of (6.3) 
and from the asymptotic solution (dashed line) with 
arg(K) = 0, arg(<|>) = —n/4. The labels (0 and (ii) on the 
graphs arefor a = 1 and 5 respectively, (b) Imaginary 
part of A'lA against Z, other parameters are as in 
Figures 6(a). (c) Real part of A'lA againstZ, other 
parameters are as in Figure 6(a). 
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